Bohr Hamiltonian with Davidson potential for triaxial nuclei 
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A solution of the Bohr Hamiltonian appropriate for triaxial shapes, involving a Davidson potential 
in P and a steep harmonic oscillator in 7, centered around 7 = it/6, is developed. Analytical 
expressions for spectra and B(E2) transition rates ranging from a triaxial vibrator to the rigid triaxial 
rotator are obtained and compared to experiment. Using a variational procedure it is pointed out 
that the Z(5) solution, in which an infinite square well potential in j3 is used, corresponds to the 
critical point of the shape phase transition from a triaxial vibrator to the rigid triaxial rotator. 
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I. INTRODUCTION 

The advent of critical point symmetries [11-3 j ma-ni- 
fested experimentally [5'-'7t| in nuclei on or near the point 
of shapephase transitions, revived interest in special so- 
lutions la, l9f of the Bohr Hamiltonian [13. While shape 
jhase transitions in nuclei have been originally found 
^ llj in the framework of the Interacting Boson Model 
12| . the first examples of critical point symmetries, the 
E(5) symmetry [l| [corresponding to the second order 
phase transition between spherical and 7-unstable (soft 
with respect to axial asymmetry) nuclei] and the X(5) 
symmetry (appropriate for the first order phase transi- 
tion between spherical and prolate deformed nuclei) , have 
been developed as special solutions of the Bohr Hamil- 
tonian, using an infinite square well potential in the 13 
degree of freedom (related to the magnitude of the de- 
formation). In E(5) the potential is independent of the 
7 degree of freedom , related to the shape of the nu- 
cleus, while in X(5) the potential is separable into two 
terms, u{f3) + v('j), the latter being a steep harmonic os- 
cillator centered around 7 = 0, corresponding to prolate 
deformed nuclei. The Z(5) solution [IJ] developed later, 
formally resembles the X(5) case in using a separable 
potential and an infinite square well potential in /3,but 
it differs drastically in using a steep harmonic oscilla- 
tor potential in the 7 degree of freedom centered around 
7 = 7r/6, corresponding to triaxial shapes. 

Triaxial shapes in nuclei have been considered for a 
long time, since the introduction of the rigid triaxial rotor 
[13, [1^ , despite the fact that very few candidates have 
been found experimentally [igI . [17| . In the framework 
of the IBM, triaxial shapes can occur in three different 
cases: 

i) In the IBM-1 framework, in which no distinction be- 
tween protons and neutrons is made, the inclusion of 
higher order (three-body) terms is needed [3 • 

ii) In the IBM-2 framework, in which protons and neu- 
trons are used as distinct entities, the inclusion of one- 
body and two-body terms suffices [20l - [2^ . 

iii) In the sdg-IBM framework, the presence of the g- 
boson also suffices [23| . 

Shape phase transitions involving rigid triaxial shapes 
have been studied recently in the IBM-2 framework [201 - 



|22| |. while in the sdg-IBM framework no transitions to- 
wards stable triaxial shapes have been found so far 23] . 

In the present work, the Z(5) solution is modified by re- 
placing the infinite square well potential in /3 by a David- 
son potential 



(1) 



where /3q corresponds to the position of the minimum of 
the potential. This solution is going to be called Z(5)-D. 
Similar studies already exist in the literature for both the 
E(5) [1 ill and X(5) [H, 113 cases. In addition, other 
potentials, like /J^" potentials [281430}. and the Morse [Slj 
and Kratz er [Sj] potentials have been used in the E(5) (7- 
unstable) [32] and X(5) (7 ~ 0, prolate deformed) [SJ] 
frameworks. 

In addition to providing easily comparable to experi- 
ment analytical solutions for the spectra and B(E2) tran- 
sition rates, the present study leads to an important by- 
product. Using a variational procedure applied earlier in 
the E(5) and X(5) frameworks [1^, [13], one can see that 
the Z(5) solution can be interpreted as corresponding to 
the critical point of a shape phase transition between a 
triaxial vibrator and a rigid triaxial rotator. 

In Sections 2 and 3 the /3 part and the 7 part of the 
spectrum are considered, while B(E2) transition rates are 
calculated in Section 4. Numerical results, including re- 
sults of the above-mentioned variational procedure, are 
shown in Section 5, while in Section 6 a brief compari- 
son to experiment is attempted. Finally, conclusions and 
plans for further work are found in Section 7. 



II. THE /?-PART OF THE SPECTRUM 



The original Bohr Hamiltonian [Toj is 
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where (3 and 7 are the usual collective coordinates, while 
Qk (fc — 1, 2, 3) are the components of angular momen- 
tum in the intrinsic frame, and B is the mass parameter. 

In the case in which the potential has a minimum 
around 7 — 7r/6, the term involving the components 
of the angular momentum can be written [13] in the 
form A{Q\ + Q\ + Q\) — 3Q^. Using this result in the 
Schrodinger equation corresponding to the Hamiltonian 
of Eq. ([2]), introducing [3] reduced energies e = 2BE/h^ 
and reduced potentials u — 2BV/h^, and assuming 
that the reduced potential can be separated into two 
terms, one depending on /3 and the other depending on 
7, i.e. u(/3,7) = u(/3) +v{'y), the Schrodinger equation 
can be approximately separated into two equations 
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where /3o is the position of the minimum of the potential. 
When plugging the Davidson potential in Eq. ([3]), the 
/3o//3^ term of the potential is combined with the [4L(L-|- 
1) — 3a^]/4/3^ term appearing there and the equation is 
solved exactly [25„ ^] , the eigenfunctions being Laguerre 
polynomials of the form 
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where T{z) stands for the F-function, n is the usual os- 
cillator quantum number (which should be distinguished 
from the wobbling quantum number n^), denotes 
the Laguerre polynomials [11], and 
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where L is the angular momentum quantum number, a 
is the projection of the angular momentum on the body- 
fixed i'-axis {a has to be an even integer [l^l), (/3^) is 
the average of (3^ over and e = ep + e^. It should be 
noticed that the separation of variables is approximate, 
since in Eq. ^ the quantity (/3^) appears, which de- 
pends on the quantum numbers L and a, appearing in 
Eq. ([3]). Therefore an approximate separation of vari- 
ables is achieved in the adiabatic limit, as in Ref. [35| . 
As a consequence, the energy relation e = + e-y is also 
approximate. 

The total wave function should have the form 
*(/3,7,^z) = a,a(/3)ry(7)2?if,a(^0, where 9, (i = 1, 2, 
3) are the Euler angles, denote Wigner functions 

of them, L are the eigenvalues of angular momentum, 
while M and a are the eigenvalues of the projections of 
angular momentum on the laboratory fixed z-axis and 
the body-fixed x'-axis respectively. 

Instead of the projection a of the angular momentum 
on the x'-axis, it is customary to introduce the wobbling 
quantum number (s^ . [36| = L — a, which labels a 
series of bands with L = n^, + 2, 4- 4, . . . (with 
riw > 0) next to the ground state band (with n^, = 0) 

M- 

Eq. ^ has been solved in the case in which u{l3) 
is an infinite well potential, the corresponding solution 
called Z(5) [l^. The spectrum is given by roots of Bessel 
functions, for which the notation has been kept the same 
as in Ref. 0, namely Es.n^,L: the ground state band 
corresponding to s = 1, = 0. 

Eq. ^ is exactly soluble also in the case in which the 
potential has the form of a Davidson potential [l^l 
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The energy eigenvalues are then (in Tiu) — 1 units) 



E^:X^^2n + a + l 
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where n = 0,1,2,. . . One can see that a formal correspon- 
dence between the energy levels of the Z(5) model and 
the present model, to which we shall refer as the Z(5)-D 
model can be established through the relation n = s — 1, 
which expresses just a formal one-to-one correspondence 
between the states in the two spectra, while the origin 
of the two quantum numbers is different, s labeling the 
order of a zero of a Bessel function and n labeling the 
number of zeros of a Laguerre polynomial. For the energy 
states the notation Es^n„,L = En+i,n„,L will be used, as 
in Refs. 0, [l^. Therefore the ground state band corre- 
sponds to s = 1 (n = 0) and riw = 0. 

In the limit /3o ^ 00 one can expand the square root 
in Eq. (|8]) and keep only the lowest order term, thus 
obtaining 



E 



A[L{L + 4) + 3n^(2L - n^)]. 
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where A is a constant, which is the spectrum of the tri- 
axial rotator obtained in Ref. (33 |. 

In the special case /3o = 0, i.e., in the case that a har- 
monic oscillator is used, one obtains a parameter-free (up 
to overall scale factors) exactly soluble model to which 
we shall refer as the Z{5)-P'^ model, in analogy to the 
E(5)-/32" [H, ^ and X(5)-/32" [3^ models. This model 
represents a triaxial vibrator. 
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III. THE 7-PART OF THE SPECTRUM 

The 7-part of the spectrum is obtained from Eq. (|4]), 
as described in Ref. [l^, by putting in it a harmonic 
osciUator potential having a minimum at 7 = 7r/6, i.e. 
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In the case of 7 ~ 7r/6 a simple harmonic oscillator equa- 
tion in the variable 7 occurs. Similar potentials and so- 
lutions in the 7- variable have been considered in [lol [39j . 

The total energy in the case of the Z(5)-D model is 
then 

E{n, Tiw, L, n;y, Pq) ^ + A[2n + I 



L{L + 4) + 3n„(2L - n„) + 9 



+ Bn^, (11) 



where n;y is the number of oscillator quanta in the 7 
degree of freedom, and Eq, A, B are free parameters. 

It should be noticed that in Eq. there is a latent 
dependence on s, L, and n^, "hidden" in the (/3^) term. 
The approximate separation of the /3 and 7 variables is 
achieved by considering an adiabatic limit, as in the X(5) 
case [i,[3f . 



IV. B(E2) TRANSITION RATES 

The quadrupole operator is given by 
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where t is a scale factor, while in the Wigner functions, 
'D^'^\ the quantum number a appears next to /i, and 
the quantity 7 — 27r/3 in the trigonometric functions is 
obtained from 7 — 27rfc/3 for fc = 1, since in the present 
case the projection a along the body-fixed x'-axis is used. 
The symmetrized wave function for Z(5)-D reads 



2L + 1 



where the normalization factor occurs from the standard 
integrals involving two Wigner functions [i^l and is the 
same as in [s^]. ce has to be an even integer [s^], while 
for a = it is clear that only even values of L are al- 
lowed, since the symmetrized wave function is vanishing 
otherwise. 



The calculation of B(E2)s proceeds as in Ref. ^T§\ and 
need not be repeated here. In the calculation of matrix 
elements the integral over 7 leads to unity [because of the 
normalization of 77(7), and taking into account that 7 in 
Eq. (|12p is fixed to the 7r/6 value, because of the steep 
potential used in 7], while the integral over /3 takes the 
form 

Il3{nt,Li, ai,nf,Lf, aj) = 



(14) 



where the (3 factor comes from Eq. (fT2l) . and the f3* 
factor comes from the volume element It is worth 

reminding, though, that a Aa = ±2 selection rule occurs, 
which results in vanishing quadrupole moments. 



V. NUMERICAL RESULTS 
A. Spectra 

The lowest bands for the Z(5)-D model are shown in 
Table 1 for the limiting parameter values /3o = (the 
Z(5)-/3^ model) and /3 — > 00 (the triaxial rotor model 
pM |39[), as well as for the intermediate value fi^ — 2 
(for illustative purposes). The levels of Z(5) [l^l are also 
shown for comparison. The bands shown are 

i) The ground state band (gsb), with [s — 1, = 0). 

ii) The quasi-71 band, composed by the even L levels 
with (s = 1,71,^, = 2) and the odd L levels with (s = 
l,n„ = 1). 

iii) The quasi-72 band, composed by the even L levels 
with (s = l,ntu = 4) and the odd L levels with (s = 
1, = 3). 

iv) The quasi-/?! band, with (s — 2,11^ — 0). 

v) The quasi- /32 band, with (s — 3, riw ~ 0). 

Since the last two bands go to infinity for /3o — > 00, the 
energy levels for /3o = 3 have been shown instead. 

In all cases B = has been used in Eq. (fTT|) . i.e., the 
term involving has been ignored. 

For all bands a uniform raising of the energies from 
the triaxial vibrator (/3o = 0) values to the triaxial rigid 
rotator (/3o — > 00) values is observed. 

A quantity being very sensitive to structural changes 
(since it is a discrete derivative of energies) is the odd- 
even staggering in gamma bands, described by quantity 



^ ' E{2+) ^ ' 

which measures the displacement of the (J — 1)+ level 
relative to the average of its neighbors, J+ and (J — 2)+, 
normalized to the energy of the first excited state of the 
ground state band, 2^ . 
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It is known that 7-soft shapes exhibit staggering 
with negative values at even L and positive values at odd 
L, while triaxial 7-rigid shapes exhibit the opposite be- 
havior, i.e., positive values at even L and negative values 
at odd L. In Table 1 it is clear that the present mod- 
els exhibit strong staggering of the triaxial type, with the 
even-L levels growing much faster with L than the odd-L 
levels. 



B. Variational procedure 

A variational procedure appropriate for locating the 
behaviour of various physical quantities at a critical point 
has been introduced ^9, H^j and applied for recovering 
the E(5) [ij and X(5) Q ground state bands from David- 
son potentials in the relevant frameworks. The method 
is applicable in cases in which one has a one-parameter 
potential spanning the region between two limiting sym- 
metries. The method is based on the fact that if a 
shape/phase transition occurs between these two sym- 
metries, the rate of change of various physical quantities 
should become maximum at the critical point [41|. The 
parameter value corresponding to the maximum, /3o,m, 
is determined for each value of angular momentum sep- 
arately. 

The variational procedure used here resembles the 
standard Ritz variational procedure of quantum mechan- 
ics (4^ . in which a trial wave function containing a free 
parameter is used, while here a potential containing a free 
parameter is used, a difference being that in the Ritz ap- 
proach the parameter is determined by minimizing the 
energy, while here the parameter is found be maximizing 
the rate of change of the relevant physical quantity. L- 
dependent potentials, like the ones occuring here, have 
been used in nuclear physics in optical model potentials 
[isl - lisj . as well as in the study of quasimolecular reso- 
nances fi^ . The method is also analogous to the variable 
moment of inertia model (VMI) 47] , in which the energy 
is minimized with respect to the moment of inertia (which 
depends on the angular momentum) separately for each 
value of the angular momentm L. 

In the present case, as seen in subsec. 2.1, the David- 
son potentials of Eq. ([T]) lead to a triaxial vibrator Z(5)- 
for (3o = 0, while they give the rigid triaxial rotator 
[13, m, [131 for /3o — >■ 00. Applying the variational proce- 
dure to the energy ratios E{L)/E{2) of the ground state 
band (s = 1, = 0) of the Z(5)-D model, where /3o is 
the free parameter serving to span the region between the 
two limiting cases, we are led to the results shown in Ta- 
ble 1, where for each value of the angular momentum L 
the location of the maximum, /3o,m, and the correspond- 
ing energy (normalized to the energy of the first excited 
state) are given. It is clear that the band determined 
through the variational procedure agrees very well with 
the ground state band of the Z(5) model. The agree- 
ment remains equally good for the s = 1, n^, = 1, 2, 3, 
4 bands, also shown in Table 1, thus indicating that the 



Z(5) model is possibly related to a shape/phase transi- 
tion from a triaxial vibrator to the rigid triaxial rotator. 



C. B(E2) transition rates 

Both intraband and interband B(E2) transition rates 
for the same models are reported in Table 2. In addi- 
tion, results for the 0(6) limit of the Interacting Boson 
Model [1^ are shown for comparison, derived from the 
expressions given in Ref. 12], the final results reading 



RgJL + 2^L) 



B{E2-{L + 2)g^ Lg) 



R 



B{E2-2g^%) 
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2 (lT5) 4:N{N + 4) ' 

B{E2;L^ Lg) 



(16) 
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^ ^ + 1) - B{E2;2g^0g) 
5(L - 1){2L + 3) {2N -L- 1){2N + L + 9) 



(17) 



L{L + 6){2L + 1) 



R. 



^_(L + 2^i) = 
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BiE2;2g^0g) 
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2(L + 7)(2L-f 3) 
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m{N + 4) 



B{E2-{L + 2)^ 
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(19) 
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(20) 



R 



AL^L-1) = 



B{E2;L^ (i- 1)7) 
B{E2;2g^0g) 



30(L + 2) {2N-L-l){2N + L + 9) 



{L - 1){L + 6){2L + 1) 



AN{N + 4) 



(21) 

In all of the above equations, N stands for the boson 
number. Numerical results for N ^ 00 are reported in 
Table 2. We remark that the 0(6) predictions for iV — )• 00 
are very similar to the ones of the rigid triaxial rotator 
[13, [il], i.e. to these of the Z(5)-D model for /3o -)■ 00. 



5 



VI. COMPARISON TO EXPERIMENT 

As seen from Table 1, one should look for nuclei having 
ground state bands charactrized by Ri/2 = E{4)/E{2) 
ratios between 2.150 and 2.667, while the 71 bandhead 
(normalized to the 2^ state) should be between 1.734 and 
2.000, the Pi bandhead (normalized in the same way) be- 
ing above 2.528 . The Xe isotopes ^^^"^■^^Xe, lying below 
the N=82 shell closure, nearly fulfil these conditions. Re- 
sults of one-parameter (/3o) rms fits are shown in Table 3, 
with a being the quality measure 



i Y.l=i{E,{exp) ~ E,{th)y 

{n-l)E{2tf 



(22) 



The overall agreement is good, with the notable excep- 
tion of the even-L levels of the quasi-71 band, which grow 
too fast, as already remarked at the end of subsec. 5.1 . 
As a result, the theoretical predictions exhibit strong tri- 
axial odd-even staggering, which is not seen experimen- 
tally. Indeed, the Xe isotopes are known [l3| to exhibit 
staggering of the 7-soft type, in contrast to the strong tri- 
axial 7-rigid staggering shown here by the theoretical val- 
ues. The only nuclei found in the extended recent search 
of Ref. [l3l to possess 71 bands with triaxial shapes are 
ii2Ru, i70Er, 1920s, i92pt^ and 232Th, ah of them located 
in the nuclear chart far from the Xe isotopes considered 
here. 

In Table 4 the existing B(E2) transition rates of the 
same nuclei are compared to the Z(5)-D model predic- 
tions for the parameter values obtained from fitting the 
spectra. No fitting of the B(E2) values has been per- 
formed. The theoretical predictions are in general higher 
than the experimental values, but in most cases lie within 
the experimental error bars, or quite near them. 



VII. CONCLUSIONS 

Z(5) [l^ is a solution of the Bohr Hamiltonian similar 
to the X(5) ^ solution, with the notable difference that 
it regards triaxial shapes (7 7r/6) instead of prolate 
deformed shapes (7 ~ 0). Predictions for spectra and 
B(E2) transition rates are parameter independent (up to 
overall scale factors). 

In the present Z(5)-D solution, the infinite square well 
m(/3) potential, used in Z(5), is replaced by the Davidson 
potential [IJ], involving a free parameter, Pq. As a result. 



Z(5)-D can cover the region between a triaxial vibrator 
and the rigid triaxial rotator [l3, [l^. In addition to 
providing easily comparable to experiment analytical so- 
lutions for spectra and B(E2) values within this wide re- 
gion, the present solution has an interesting by-product. 
Using a variational procedure [1^ ^\ it is pointed out 
that the Z(5) solution corresponds to the critical point of 
the shape phase transition from a triaxial vibrator to the 
rigid triaxial rotator. However, the Z(5) solution is not a 
special case of Z(5)-D, obtained for a specific parameter 
value, or a limiting case of Z(5)-D. Using the Davidson 
potential one can cover the whole way from triaxial vi- 
brator to triaxial rotator, but one cannot get the critical 
point as a special case. This is due to the shape of the 
Davidson potential, which is not fiat, as the potential is 
expected to be at the critical point. The same situation 
has occured in the ESD model [5l|, in which the David- 
son potential is used in order to interpolate between a 
vibrator and the prolate axial rotator with 7 « 0. Using 
the ESD model one can obtain very good fits of many 
nuclei from the prolate rotator limit down to close to the 
critical point, but one cannot describe the nuclei very 
close to the critical point [Slj . 

Concerning the separation of variables which allowed 
for analytical solutions, a potential of the form u(/3) -I- 
^(7) has been used, bringing in the approximations used 
in X(5) '2^. These approximations can be avoided in two 
ways: 

i) Using potentials of the form u{f3) -\- v{'y)/f3'^, which 
are known [8| to allow for exact separation of variables 
without any approximations. 

ii) Using the powerful techniques of the Algebraic Collec- 
tive Model (53 - t54| . which allow for the exact numerical 
diagonalization of any Bohr Hamiltonian. 

The first path has been used for a detailed study of the 
Davidson potential plugged in the Bohr Hamiltonian for 
7 sa 51]. The main advantage of this solution is that 
all bands are treated on equal footing with respect to the 
influence of the 17(7) potential, while in the present solu- 
tion only the quasi-7 bands are affected. A similar study 
for 7 7r/6 case would be interesting. The analytical 
solution and a brief comparison to experiment in the Os 
region has already been given in Ref. (ssj . 

The second path has been recently used for the de- 
scription of a triaxial symmetry top [54], as well as for 
the study the onset of rigid triaxial deformation [5^ . Fur- 
ther investigations of triaxial shapes using this powerful 
tool should also be revealing. 
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TABLE I: Energy spectra of the Z(5)-D model (for different values of the parameter /3o), and for the Z(5) model [T3|- l3o — 
corresponds to the Z(5)-/3^ model (a triaxial vibrator), while /? — i- oo is the rigid triaxial rotator fl3.[T5l]. The notation I/s,n„ is 
used. All levels are measured from the ground state, Oi,o, and are normalized to the first excited state, 2i,o- See subsec. 5.1 for 
further discussion. In addition, the energy levels resulting from the variational procedure of subsec. 5.2 are reported (labelled 
by "var"), along with the parameter values /3o,m at which they are obtained. See subsec. 5.2 for further discussion. 



l3o 





2 


oo 








/3o 





2 


oo 


















var 


Z(5) 














var 


Z(5) 


Oi.o 


0.000 


0.000 


0.000 






0.000 


















2i,o 


1.000 


1.000 


1.000 






1.000 


















4i.o 


2.150 


2.521 


2.667 


1.375 


2.341 


2.350 


















6i,o 


3.353 


4.424 


5.000 


1.474 


3.956 


3.984 


















8i,o 


4.579 


6.596 


8.000 


1.562 


5.819 


5.877 


















10i,o 


5.817 


8.957 


11.667 


1.640 


7.915 


8.019 


















12i,o 


7.063 


11.450 


16.000 


1.713 


10.237 


10.403 


















14i,o 


8.313 


14.039 


21.000 


1.780 


12.781 


13.024 


















16i,o 


9.566 


16.698 


26.667 


1.843 


15.544 


15.878 


















18i,o 


10.821 


19.410 


33.000 


1.902 


18.523 


18.964 


















20i,o 


12.077 


22.163 


40.000 


1.960 


21.719 


22.279 
































2i 


2 


1.734 


1.932 


2.000 


1.336 


1.833 


1.837 


3i,i 


2.343 


2.807 


3.000 


1.392 


2.586 


2.597 


4i 


2 


3.649 


4.930 


5.667 


1.496 


4.386 


4.420 


5i,i 


3.791 


5.177 


6.000 


1.507 


4.597 


4.634 


6i 


2 


5.281 


7.917 


10.000 


1.607 


6.981 


7.063 


7i,i 


5.169 


7.703 


9.667 


1.600 


6.790 


6.869 


8i 


2 


6.791 


10.898 


15.000 


1.697 


9.713 


9.864 


9i.i 


6.511 


10.333 


14.000 


1.681 


9.182 


9.318 


lOi 


2 


8.234 


13.874 


20.667 


1.776 


12.615 


12.852 


lli.i 


7.832 


13.035 


19.000 


1.754 


11.778 


11.989 


12i 


2 


9.635 


16.847 


27.000 


1.846 


15.703 


16.043 


13i,i 


9.140 


15.788 


24.667 


1.822 


14.581 


14.882 


14i 


2 


11.008 


19.818 


34.000 


1.911 


18.986 


19.443 


15i,i 


10.438 


18.579 


31.000 


1.885 


17.593 


18.000 


16i 


2 


12.360 


22.787 


41.667 


1.972 


22.468 


23.056 


17i,i 


11.730 


21.399 


38.000 


1.944 


20.815 


21.341 


18i 


2 


13.698 


25.755 


50.000 


2.030 


26.154 


26.884 


19i,i 


13.017 


24.241 


45.667 


2.001 


24.248 


24.905 


20i 


2 


15.024 


28.721 


59.000 


2.085 


30.046 


30.928 



5l,3 
7l.3 
9l,3 
lll,3 
13l,3 
15l,3 
17l,3 
19l,3 



4.939 
6.699 
8.313 
9.841 
11.314 
12.747 
14.152 
15.536 



7.268 
10.711 
14.039 
17.289 
20.486 
23.642 
26.768 
29.871 



9.000 
14.667 
21.000 
28.000 
35.667 
44.000 
53.000 
62.667 



1.585 
1.692 
1.780 
1.856 
1.925 
1.989 
2.049 
2.105 



6.406 
9.537 
12.781 
16.180 
19.752 
23.509 
27.460 
31.611 



6.476 
9.683 
13.024 
16.536 
20.237 
24.137 
28.241 
32.553 



4l,4 
6l,4 
8l,4 
10l,4 
12l,4 
14l,4 
16l,4 
18l,4 

20i,4 



4.066 
6.221 
8.075 
9.773 
11.374 
12.910 
14.399 
15.853 
17.281 



5.663 
9.753 
13.541 
17.143 
20.618 
24.003 
27.320 
30.586 
33.811 



6.667 
13.000 
20.000 
27.667 
36.000 
45.000 
54.667 
65.000 
76.000 



1.526 
1.665 
1.767 
1.853 
1.928 
1.996 
2.059 
2.118 
2.174 



5.012 
8.644 
12.282 
16.021 
19.904 
23.953 
28.182 
32.601 
37.217 



5.056 
8.767 
12.508 
16.372 
20.396 
24.598 
28.991 
33.581 
38.373 



/3o 







/3o 













/3o,m 


var Z(5) 












var Z(5) 


02,0 


2.528 


5.921 


12.274 


3.913 


03,0 


5.055 


11.842 


24 


548 


9.782 


22,0 


3.528 


6.921 


13.274 


5.697 


23,0 


6.055 


12.842 


25 


548 


12.343 


42,0 


4.678 


8.442 


14.903 


7.962 


43,0 


7.205 


14.363 


27 


177 


15.506 


62,0 


5.881 


10.345 


17.110 


10.567 


63,0 


8.408 


16.266 


29 


384 


19.059 


82,0 


7.107 


12.517 


19.835 


13.469 


83,0 


9.634 


18.439 


32 


109 


22.933 


102,0 


8.345 


14.878 


23.015 


16.646 


103,0 


10.873 


20.799 


35 


289 


27.103 


122,0 


9.590 


17.371 


26.588 


20.088 


123,0 


12.118 


23.292 


38 


862 


31.552 


142,0 


10.840 


19.960 


30.497 


23.788 


143,0 


13.368 


25.881 


42 


771 


36.272 


162,0 


12.093 


22.619 


34.692 


27.740 


163,0 


14.621 


28.540 


46 


966 


41.258 


182,0 


13.348 


25.331 


39.129 


31.942 


183,0 


15.876 


31.253 


51 


403 


46.504 


202,0 


14.605 


28.084 


43.772 


36.390 


203,0 


17.132 


34.005 


56 


046 


52.007 
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TABLE II: Intraband and interband B(E2) transition rates, normalized to the one between the two lowest states, B(E2;2i,o — >■ 
Oi,o), are given for Z(5)-D model (for different values of the parameter j3o), for the Z(5) model fl3|, and for the 0(6) limit of 
the Interacting Boson Model [l^. /3o = corresponds to the Z(5)-/?^ model (a triaxial vibrator), while /3 ^ oo is the rigid 
triaxial rotator ITsI]. The notation Ls,n„ is used, while the initial state is labelled by (i), and the final state by (/). 



/3o 







2 




00 






A) 







2 




00 




rW 


r(/) 










0(6) 


Z(5) 


/-{») 


r(/) 










0(6) Z(5) 


2l,0 


Oi,o 


1.000 


1.000 


1 


000 


1.000 


1.000 


22,0 


02,0 


1.480 


1.284 


1 


000 


0.774 


4i,o 


2i,o 


1.834 


1.493 


1 


389 


1.429 


1.590 


42,0 


22,0 


2.440 


1.865 


1 


389 


1.192 


6i,o 


4i,o 


2.919 


2.041 


1 


731 


1.667 


2.203 


62,0 


42,0 


3.647 


2.477 


1 


731 


1.643 


8i,o 


6i,o 


3.955 


2.497 


1 


912 


1.818 


2.635 


82,0 


62,0 


4.746 


2.955 


1 


912 


1.975 


lOi.o 


8i,o 


4.976 


2.934 


2 


024 


1.923 


2.967 


102,0 


82,0 


5.804 


3.398 


2 


024 


2.242 


12i,o 


lOi.o 


5.989 


3.370 


2 


100 


2.000 


3.234 


122,0 


102,0 


6.844 


3.836 


2 


100 


2.466 


14i,o 


12i,o 


6.999 


3.810 


2 


155 


2.059 


3.455 


142,0 


122,0 


7.874 


4.277 


2 


155 


2.660 


16i,o 


14i,o 


8.007 


4.257 


2 


197 


2.105 


3.642 


162,0 


142,0 


8.897 


4.723 


2 


197 


2.829 


18i,o 


I61.0 


9.014 


4.708 


2 


230 


2.143 


3.803 


182,0 


162,0 


9.915 


5.175 


2 


230 


2.980 


20i,o 


18i,o 


10.020 


5.164 


2 


256 


2.174 


3.944 


202,0 


182,0 


10.931 


5.630 


2 


256 


3.115 



4l,2 

61,2 
81,2 

10l,2 
12l,2 
14l,2 
16l,2 
18l,2 

20i,2 



2l,2 
4l,2 

61,2 
81,2 

10l,2 
12l.2 
14l,2 
16l,2 
18l,2 



0.912 
1.583 
2.816 
4.038 
5.231 
6.395 
7.535 
8.655 
9.759 



0.682 
0.986 
1.617 
2.215 
2.786 
3.339 
3.877 
4.406 
4.927 



0.595 
0.734 
1.051 
1.278 
1.446 
1.574 
1.674 
1.755 
1.822 



0.873 
1.240 
1.462 
1.614 
1.726 
1.813 
1.882 
1.938 
1.985 



0.736 
1.031 
1.590 
2.035 
2.394 
2.690 
2.938 
3.151 
3.335 



3i,i 


2l,2 


5i,i 


4l,2 


7i,i 


61,2 


9i,i 


81,2 


lli,i 


10l,2 


13i,i 


12l,2 


15i,i 


14l,2 


17i,i 


16l,2 


19i,i 


18l,2 


3i,i 


4i,o 


5i,i 


61,0 


7i,i 


81.0 


9i,i 


lOi.o 


lli,i 


12i,o 


13i,i 


14i,o 


15i,i 


16i,o 


17i,i 


18i,o 


19i,i 


20i,o 



2.731 
1.991 
2.183 
2.247 
2.265 
2.265 
2.258 
2.247 
2.236 



2.006 
1.244 
1.262 
1.241 
1.214 
1.189 
1.168 
1.149 
1.133 



1.786 
0.955 
0.851 
0.746 
0.657 
0.585 
0.526 
0.478 
0.437 



1.190 
0.434 
0.231 
0.145 
0.100 
0.073 
0.056 
0.044 
0.036 



2.171 
1.313 
1.260 
1.164 
1.069 
0.984 
0.910 
0.846 
0.790 



2l,2 
4l,2 

61,2 
81,2 

101.2 

12l,2 
14l,2 
16l,2 
18l,2 

20i,2 



2i,o 
4i,o 
61,0 
81,0 
10i,o 
12i,o 
14i,o 
16i,o 
18i,o 
2O1.0 



1.865 
0.459 
0.292 
0.211 
0.165 
0.135 
0.114 
0.099 
0.087 
0.078 



1.520 
0.323 
0.187 
0.127 
0.094 
0.074 
0.061 
0.052 
0.045 
0.040 



1.429 
0.273 
0.143 
0.088 
0.059 
0.043 
0.032 
0.025 
0.020 
0.017 



1.429 
0.794 
0.579 
0.462 
0.386 
0.332 
0.292 
0.261 
0.236 
0.215 



1.620 
0.348 
0.198 
0.129 
0.092 
0.069 
0.054 
0.043 
0.035 
0.030 



1.618 
1.449 
1.351 
1.287 
1.243 
1.211 
1.186 
1.167 
1.151 



1.147 
0.917 
0.796 
0.724 
0.676 
0.643 
0.619 
0.601 
0.587 



1.000 
0.714 
0.556 
0.455 
0.385 
0.333 
0.294 
0.263 
0.238 



0.476 
0.430 
0.374 
0.327 
0.291 
0.261 
0.237 
0.216 
0.199 



1.243 
0.972 
0.808 
0.696 
0.614 
0.551 
0.507 
0.459 
0.425 



5i,i 
7i,i 
9i,i 
lli.i 
13i,i 
15i,i 
17i,i 
19i,i 



3i.i 
5i,i 
7i,i 
9i,i 
lli,i 
13i,i 
15i,i 
17i,i 



1.667 
2.891 
4.061 
5.191 
6.292 
7.373 
8.440 
9.496 



1.147 
1.778 
2.338 
2.865 
3.374 
3.873 
4.366 
4.856 



0.955 
1.310 
1.535 
1.690 
1.802 
1.887 
1.954 
2.007 



0.955 
1.319 
1.528 
1.668 
1.771 
1.850 
1.913 
1.965 



1.235 
1.851 
2.308 
2.665 
2.952 
3.190 
3.392 
3.566 



9 



TABLE III: Comparison of theoretical predictions of the Z(5)- 
D model, labelled by the relevant /3o value, to experimental 
spectra of ^^*Xe [3, ""Xe [11], and "^Xe In each col- 
umn all energies are normalized to the energy of the relevant 
2f state. The quality measure a of Eq. (|22p is used. See Sec. 
6 for further discussion. 





^^«Xe 


^^«Xe 


^^•^Xe 


^^"Xe 


^^^Xe 


^^^Xe 




exp 


Po = 1.32 


exp 


/3o = 1.11 


exp 


/3o = 


4i,o 


2.333 


2.323 


2.247 


2.255 


2.157 


2.150 


Dl,0 




O.oUO 


o.OZ t 


O.DZl 


o. IDO 


o.oOo 


8i,o 


5.674 


5.372 


5.031 


5.040 






lOi.o 


7.597 


6.986 


6.457 


6.489 






12i,o 






7.867 


7.956 






14i,o 






9.458 


9.434 






2l,2 


2.189 


1.830 


2.093 


1.793 


1.944 


1.734 


4l,2 


3.620 


4.180 


3.373 


3.961 


2.940 


3.649 


6l,2 


5.150 


6.284 










3i,i 


3.228 


2.555 


3.045 


2.471 


2.701 


2.343 


5i,i 


4.508 


4.360 


4.051 


4.125 


3.246 


3.791 


7i,i 


6.165 


6.138 










02,0 


3.574 


3.452 


3.346 


3.028 


2.771 


2.528 


22,0 


4.515 


4.452 










a 




0.495 




0.297 




0.422 



TABLE IV: Comparison of theoretical predictions of the Z(5)- 
D model, labelled by the relevant /3o value, to experimental 
B(E2) values of ^^**Xe :4gl and "^Xe [Ml- In each column all 
B(E2)s are normalized to the 2^ — >■ 0^ transition. See Sec. 6 
for further discussion. 







i^«Xc 


^^«Xc 


i^^Xe 


^a^Xe 


r{») 


r(/) 


exp 


Po = 1.32 


exp 


/?o = 


4i,o 


2i,o 


1.468 ±0.201 


1.648 


1.238 ±0.180 


1.834 


6i,o 


4i,o 


1.940 ± 0.275 


2.464 






8i,o 


6i,o 


2.388 ± 0.398 


3.228 






2l,2 


2i,o 


1.194 ±0.187 


1.673 


1.775 ±0.288 


1.865 


2l,2 


Oi.o 


0.016 ±0.002 


0.000 


0.003 ± 0.001 


0.000 



